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Abstract. For a general quantum theory that is describable by a path integral 
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whose outcomes give predictions that are nearly identical to the given quantum theory. 
The model is neither local nor causal in spacetime, but is both local and causal is in 
a non-observable path space. The probabilistic nature of the squared wavefunction 
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detectable via experiment. Finally, we discuss the physical implications of the model. 

Keywords: quantum theory, quantum mechanics. Born rule, signal processing, 
threshold process, path integral, causality 


PACS numbers: 03.65.Ta 



An Accumulative Model for Quantum Theories 

1. Introduction 


2 


The paradoxical, apparently indeterministic nature of quantum theory is common 
knowledge, and numerous attempts have been made to provide a deterministic, causal 
basis for the theory. One possible approach is to admit the possibility of causes outside 
of space and time. Bohm takes this point of view, and identihes spacetime events as 
’’unfoldings” of a more fundamental ’’implicate order” that is manifested within ’’pre¬ 
space” [1]. Frescura and Hiley, building on Bohm’s idea, have developed an algebraic 
representation of pre-space dynamics [2]. A somewhat different tack is taken in [3], which 
essentially proposes a conceptual model of pre-space and a deterministic, local (but 
statistically random) dynamics within that pre-space that, when “unfolded”, produces 
quantum particle transmission. The model derives from an analogy to signal detection 
in wireless communications: particle detection is represented as the outcome of a signal 
accumulation process which occurs in spacetime augmented by an extra, non-spacetime 
dimension (referred to as the a-dimension). The quantum wavefunction corresponds 
to in-phase and quadrature-phase components of an amplitude and phase-modulated 
carrier signal field that is present throughout spacetime augmented by the a-dimension. 
The location of particle detection is determined when an accumulated signal reaches a 
threshold (so that attaining the threshold effects the “unfolding”). The paper proves 
that the Born probability rule is a mathematical consequence: however, the paper gives 
no explanation of the origin or formation of the carrier signal field required for the 
model. 

The current paper provides a more comprehensive interpretation of quantum 
probabilities than [3] by taking a related, but somewhat different approach. The 
approach is based on the observation that both quantum mechanics and quantum 
field theory may be derived from a path integral formalism. We conjecture that 
path integrals correspond to a universal physical process which essentially performs a 
numerical integration. As in the previous paper, this process unfolds in a non-spacetime 
dimension, and the observable universe is the outcome of the process upon attaining a 
threshold. 

The paper is organized as follows. Section presents a simplihed preliminary 
mathematical model which illustrates the basic model structure. We demonstrates 
the model’s ability to generate quantum probabilities both theoretically and with 
simulations. Section gives a more detailed model which is designed to conform 
more closely with the hypothesized physical processes involved. Section discusses 
the possibility of experimental verihcation of the model; and Section gives a summary 
discussion. For the sake of completeness, the Matlab/Octave source code used in our 
simulations is given in Section 
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Let lA represent the space of all possible configurations of the observable universe. We 
emphasize that any u E U expresses the entire conhguration of the universe over all 
times, not just its conhguration at a single time. We do not need to specify whether we 
are employing a quantum-mechanical or held-theoretic representation of the universe’s 
conhguration space - our argument does not depend on the specihc nature of U. 

In both quantum-mechanical or the held-theoretic representations of U, the 
wavefunction can be expressed in terms of a path integral : W —)■ C of the form:: 






76r„ 


( 1 ) 


where r„ is a space of paths corresponding to the conhguration u, and is the 

action associated with the path 7. Here we have used summation notation to facilitate 
the connection with simulations that we will describe later. We shall suppose that |r„| 
is independent of u, so that |r^j| = |r|/|W| where T = U„r„. We also suppose that the 
{Tu}u&a are disjoint, which implies that for every 7 G L there exists a unique eU 
such that 7 G 

The path integral is associated with a probability distribution: 


Ps{u) ^ 




( 2 ) 


The fact that this probability is written in terms of a summation (or integral) suggests 
that some sort of accumulation process could be involved. The main purpose of 
this paper is to show that such an interpretation is indeed feasible, and provides a 
simple, plausible explanation of the hidden dynamics that give rise to quantum theories. 
Preliminarily, we note that our interpretation must address two issues: 

• Why is probability obtained from a squared complex amplitude? 

• What physically corresponds to the division in ([^? 

In the following, we give what we believe to be satisfactory answers to these two 
questions. 

We dehne an accumulation process as follows. Given the sequence of paths 71,72,... 
in T, we dehne an accumulated amplitude Ak {K E Z+) as: 


Ak = 


E 


K 

k=l 


g*>5(7fe) 


( 3 ) 


One possible interpretation of each factor is as the phasor representation | 3 ] 

of an oscillation (of unknown frequency) which depends on 7^. The summation then 
corresponds to the complex amplitude of a harmonic oscillator (with the same frequency) 
that is successively perturbed by these oscillations. 

Although we are using discrete notation, the sequence {7^,} should be thought of as 
a discrete approximation of a path-valued function of a continuous index, corresponding 
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to a continuously-varying path within the space V of all possible paths. The continuous 
index corresponds the the a-dimension introduced in [3]: and the variation within T 
corresponds to an evolutionary process within this dimension which uniformly samples 
r over the long term. Note that as 7 ^ varies, the corresponding state of the universe 
Uk = also varies. In the process we will define, the accumulated amplitude grows to 
reach a fixed threshold at a particular index K, at which point uk gives the configuration 
of the observable universe. 

In order to obtain the probabilities ([^ via this process, we impose additional 
conditions on the sequences { 7 ^,} and {uk} as follows. 

(a) There exists N ^ 1 and M ^ 1 such that UkNM+i = UkNM +2 = ■ ■ ■ = 

U{k+i)NM,^k G Z>o; 

(b) For each k e Z>o, the sequence {'JkN+i,'JkN+ 2 , ■ • •, 7(fc+i)Af} uniformly samples 

(c) The sequence {unm,U 2 nm, ...} is mixing [5] and uniformly samples U. 

These conditions correspond to a situation where { 7 ^} varies throughout T such that 
the sequence { 7 ^} uniformly samples each T^ that it visits before passing on to the next 
r„. In this simple model, the dwell time within each r„ visited is the constant N: in 
our subsequent model, this assumption will be relaxed. The significance of M will be 
explained later. 

Let rjk {k = 1,2,...) be a sequence of independent, identically distributed (i.i.d.) 
complex-valued random variables with zero mean and finite variance, and define: 

K 

V\k/N] (4) 

k=l 

Finally, given 0 > 0, we define the threshold index as the random variable: 

Kq = min(/c||^'^| < 0 and l^'^l > 0). (5) 

Given the above conditions and definitions, we have the following result: 

Proposition: As N, M, 6 —)■ 00 , we have 

In other words, the probability distribution on U at the stopping time defined by 
attaining the threshold 6N\fM agrees with the probability distribution ([^ obtained 
from the path-integral formalism. 

The proof of this proposition is similar to that given in [3]. Notice that Q can be 
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A 


K 


N 


■KN 


9N^/M 9N^/M 




J ^ 5" ('7 { A; — 1) iV+n ) 


k=l 


^n=l 


\u\ 


K 


N^oo 




(7) 

( 8 ) 


/ IK/M\ ^ M 


M 
^iri V 


'' m=l 


K 




(9) 


A' 

^ r, 


The proof is based on the fact that can be approximated in distribntion as a 

Brownian motion B{a) in C with absorbing bonndary at 71 = 1, where a = qt§m- 
For any fixed a, near the boundary the probability density of an absorbing Brownian 
motion is proportional (to first order) to the distance from the boundary. This can be 
used to show that for any K, the probability = K\uk = u) is approximately 

proportional to E\\r]K'^{u)^^], which is proportional to |T(m) 7 Since for P{uk = u) is 
independent of u when 1 ^ iF < it follows that P{Kgj^^ = iF & = u) 

is proportional to |T(n)p, and summing over K gives the desired result. 

Figure [T] shows the results of simulations of the model specified by conditions 
0-0 and equations 0-(§. The simulations were performed on a discrete system 
with 11 possible states. To shorten computational time, the simulation was based on 
equation 0 rather than performing the full computation 0 on a path-by-path basis. 
The random variables {unm,U 2 NM, ■ ■ ■} referred to in Q were generated uniformly 
randomly. The curves show the difference between the simulated probabilities and 
actual probabilities for two different probability distributions |Tp, for different values 
of the threshold 9. The errors are shown on the j/-axis, versus the actual probability 
values which are shown on the x-axis. As 9 increases, the errors decrease: for 9 = 40, 
the maximum error is under 5 percent. The pattern of error apparently depends on 
the type of probability distribution being modeled. However, in both cases the larger 
probabilities are underestimated, and there is a range of intermediate probabilities that 
are overestimated. These phenomena may possibly enable an experimental test of the 
model: this possibility is explored further in Section 


3. Refined model 

The process we have presented above has some seemingly artihcial features: 

• Why should {uk} remain constant for intervals of size MNl 

• What is the physical signihcance of the r^^’s? 
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disinbuinn of Ivt Simulalion error 



UrMform dislnbulion of Simulatron error 





Figure 1. Deviations of computed probabilities from quantum values, for simulated 
preliminary accumulation model with 6 = 10,20,30,40 and M = 10000, where {rjk} 
are i.i.d. standard normal random variables. Each simulation was run 100,000 times. 
All simulations used 11 configurations u. For the figure on the left, the \ijj{uj)\ oc 
ji 0 = 0i • • ■ 1 10); while for the figure on the right, \ip{uj)\‘^ oc j (j = 0,..., 10). 


As to the first point, instead of supposing that {uk} remains constant on intervals of 
size N, we may suppose that {uk\ varies slowly with fc, so that 

p{uk+i ^Uk) = 0 . ( 10 ) 

Supposing that { 7 A:}fc=i. 2 . 3 .... is generated by a Markov process, it is reasonable to suppose 
that residence times in each u state visited are (approximately) i.i.d. geometrical random 
variables. This is because under reasonable conditions, hitting times in Markov chains 
are asymptotically exponentially distributed [B]. (The geometrical distribution is the 
discrete analog of the exponential distribution.) Accordingly, we may modify the model 
by replacing the constant M with a geometrically-distributed random variable with the 
same mean. 

As to the second point, we must recognize that we have failed to account for the fact 
that in practice we never measure the state of the entire universe, but only a subsystem. 
So we must take into account the effect of variations in the external system during 
the accumulation process. Accordingly we let be the possible states of the measured 
subsystem, while ff' denote the possible states of the universe external to the measured 
subsystem. Thus we may represent any element u eU uniquely as m = {w,w'), where 
w E fl and w' E ff'. 

We suppose that any path in T can be factored into a part for and a part 
for more precisely, that there are path spaces C and C' respectively such that 
any 7 G T can be decomposed as 7 = (c, d) where c E C,c' E C, and such that 
We dehne = {c\wc = w}, and suppose (as in the simple model) that 
\Cw\ is independent of ta G fl, so that \Cw\ = |C|/|f2| \/w. We similarly dehnine C'^,, 
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and suppose \C!^,\ = |C'|/|fi'| Ww'. Finally, we suppose that the action S is additive: 


Silk) = 

Sick) + Sic}). 

From this it follows that we may write: 




T(n) = T((w,t(;')) = '^(w)0(t(;'). 

(11) 

where 

f)iw\ 

, ^ M ^ •#>(»') ^ IS E 

' ' cec™ 1 1 c'ec„/ 

(12) 


We may also rewrite ([^ as 


Ak = 




(13) 


We now postulate the existence of a Markov chain {(ci, c'^), (c2, C2),...} that satisfies 
the following properties. Define inductively a sequence of random times {X^} such that 

Xo = 1; Xk+i = min(j|u;' ^ 

We suppose the Markov chain has transition probabilities such that Wj 7^ 

Wj 7^ This supposition reflects the assumption that the external state varies 

more rapidly than the observed state, which is reasonable since the external state is 
much, much larger and has many more possibilities for variation. In this case, it is 
possible to define inductively a sequence of random times { 2 k} such that Zq = 1 and 
Zk+i = 7^ According to these definitions, the state w' does not change 

on each time interval [Afc, Xk+i — l], and the state w does not change on each time interval 
[Xz^i Xz^^j^-i]. We also suppose the paths vary much faster than the states, so that the 
space is uniformly sampled on the time interval [Afc, Xk+i — 1]. 

Based on the Markov chain described in the previous paragraph, we may formulate 
the following model assumptions: 

(A) There exists a A^ S> 1 and a sequence {^1, ^2, • • •} of i.i.d. geometrically-distributed 

random variables with E[^k] = Xf, such that = • • • = ViF G 

Z>o, where Xq = 0 and Xk = ■Cfc, K > 1; 

(B) There exists a M S> 1 an a sequence {Ci, C2,...} of i.i.d. geometrically-distributed 

random variables with E[(k] = Af, such that wx^^+i = wxzj^+2 = • • • = 
'^^zk+1 ^ ^^0’ where ^0 = 0 and Zk = Y.k=i Ck,K>l] 

(C) For each K e Z>o, the sequences c!xk+2^ ■ • • > ^'xK+ik^ [cxk+i^ cxii+2, ■ ■ ■ Xxk 

uniformly sample C'x^^^ and C, respectively; 

(D) The sequences c',.^+2, • • • J {cxk+iXXk+ 2 , ■ ■ ■ XxK+ik) are 

statistically independent; 

(E) The sequences {w'x^^w'x^,...} and [wxz^^wxz^, ■ ■ ■} are mixing, and uniformly 
sample D and D' respectively. 
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A' 


Zk 


K —1 ^k + l 1 ^m + 1 




0nVm 0n^ 

1^11^1 ^m+1 ,/ \rK( ' \ 


d^\C\\C'\ t, 


k=0 m=Z^. 


( 14 ) 




K-l 


e\c\\o 


- V’(wa', 


^fc+i-i 

E 


^m+l 


k=0 


^"+1 ^/M ^ N 

m=Zk 








e\c\\c' 

where the approximation holds for large N and 


(15) 


^m,fc — 





+m 


N 


w 




fc-l+^ 


(16) 


Notice the similarity between ([^ and (15). Instead of a summation over M, there 
is a summation over which has expectation M. Within this summation, instead 
of the mean-zero i.i.d. random variables {rjk}, we now have {rjk^m} given by the 


complicated expression (16). By assumption, the variables C,k/M and izu-i+m/N are 


independent, and have expectation 1; while the additional complex factor 0 \ w';yz + 
will vary randomly with mean zero as the process evolves. If we assume that {rjk^m} 
are (approximately) i.i.d. mean-zero random variables, then (15) and are virtually 


identical, except that Cfc in (15) replaces M in ([^. However, E[C^k] = M] and 
conditioning on the different possible values of C,k-, we may obtain the same result that 
the probability density for wkq is given by \%Ij{w)\^. 

Figure shows results of simulations of the refined model specified in (|^-([^. A 
system with 31 discrete states was simulated, and the states’ probabilities were chosen 
according to the sinusoidal wavefunction shown in the picture. The transition between 
states w was determined according to a Markov chain that produced a mean dwell 
time of M, followed by a transition to one of the four nearest-neighbor states with 
equal probability 1/4. Parameters used were M = 625 and 6 = 10. The figure shows 
very close agreement between quantum-theoretic probabilities and those obtained from 
simulation. Deviations are shown in more detail in Figure]^ for different values of M 
and 9. Small IV’P’s are consistently overestimated, and large are underestimated. 

Deviations between simulation and quantum theory decrease with increasing M and 
6 *, so that the model probabilities apparently converges to quantum-theoretic values as 
M,6 ^ oo. 
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Wavefunction if* used in sinnulation Simulated versus quantum-theoretic probability 



index i 


Figure 2. (Left) Sinusoidal “wavefunction” used in simulation. 31 configurations were 
used with probabilities as shown. (Right) Simulation results compared to theory for 
0 = 10, M = 625. Computed probabilities are based on 10 million repetitions. 


Effect of threshold parameter 8 on simulation probability 



X ^0'"* accumulation length M on simulation probability 


s -c 


$ 


* * $ 


* M=25 

+ M=125 

t*t=625 


+ 

+ 


* 

* 

rfi 


+ 

* 


+ 

* 




$ 


+ +=!= 

K 


0.01 


002 


003 004 

IVl^ 


005 


006 


007 


Figure 3. Deviations of computed probabilities from quantum values for simulated 
adjusted accumulation model, for different values of the accumulation length M and 
threshold parameter 6 (as specified in the figure titles). All computed probabilities are 
based on 10 million repetitions. 


4. Proposed Experimental Test 

In the above model, quantum probabilities are generated by an accumulative process 
which essentially performs a stochastic approximation to the quantum path integrals. 
In the previous section we showed that hnite values of 6 and M introduced deviations 
from quantum-theoretical probabilities. In both cases, the deviations are positive for 
small probabilities, but negative for large probabilities. 

Another possible source of numerical error, which we did not model in the 
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Ti - <^772 + 1 


^iS{Cn) 




(17) 


which was used in (14). If we suppose there is a random error of constant variance 


in this approximation, then by carrying through the computations it can be shown 
that probabilities turn out to be proportional to rather than 

This produces a deviation from theoretical probabilities that decreases linearly with 
increasing probability density. So the deviations from quantum-theoretic probabilities 
due to this effect reinforce the deviations already discussed. 

We may conclude that numerical approximation effects should introduce a deviation 
from quantum-theoretic probabilities that for larger probabilities decreases roughly 
linearly with increasing probability density. Unfortunately, since the parameters of the 
process are not directly accessible, it is not possible to predict the size of the deviations. 


5. Discussion 

This construction provides a conceptually simple solution to many conundrums of 
quantum theory. It accounts for all quantum paradoxes, since it yields the same 
probabilities as quantum theory (to a close approximation). The model requires no 
distinction between observer and observed, because the probabilistic signihcance of the 
wavefunction is a consequence of the model, rather than an extraneous assumption that 
is made to match theory with experiment. In particular, the “measurement problem” is 
no longer a problem: what is perceived as a “collapse of the wavefunction” corresponds to 
the fact that one particular state of the universe is selected as a result of the thresholding 
process. 

According to the model, the universe is not at all what it appears. It indicates 
that causality is an illusion, and that apparent “cause and effect” relationships are 
correlations in the outcome of an inaccessible process that occurs outside of spacetime. 
The Big Bang does not account for the “origin” of the universe, because it is also part of 
the outcome of an extra-dimensional process which produces past, present, and future 
together as an entirety. (The model thus seems to imply that the universe will have hnite 
duration.) The vacuum is not a “boiling sea of virtual particles and antiparticles,” [7j as 
quantum held theories seem to imply, but only appears so because of the accumulation 
process through which the observable universe is actualized. 

We mentioned in the introduction that the “pre-space” approach proposed and 
developed by Bohm, Hiley, and others bears some similarity to our approach. However, 
our portrayal of pre-space is radically different from that envisioned by Bohm et. ah . 
Perhaps the best way to see this is to compare motivating analogies. 

The original inspiration for our model was the example of a cellular phone which 
accumulates a pilot signal broadcast by a base station as the phone is carried in the 
user’s pocket. When the signal accumulation reaches a certain threshold, a detection is 
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logged. The users location at the moment of detection is determined by the process of 
signal accumulationbut no record of his past motion may be seen in the final outcome. 
Still, the outcome reflects the process in that the detection location is more likely to be 
at a location where the signal is strong. In other words, the legacy of the process of 
signal accumulation is seen in the probability distribution of the observed outcome. 

On the other hand, Bohm in [1] describes an experiment in which a droplet of 
dye is introduced into a viscous fluid, the fluid is stirred, and the process is repeated 
several times. When the fluid is stirred in the reverse direction, the droplets reappear 
one by one. These droplets represent the unfolded order that is evidenced in spacetime 
events. Fresca and Hiley take this illustration as a jumping-off point in their portrayal 
of quantum processes in terms of successive enfolding and unfolding. Thus spacetime 
events are conceived as manifestations of an ongoing process. Clearly this is very different 
from our description of a process from which all the entire history of the universe springs 
full-blown into existence, like Athena springing from the head of Zeus. 

At this point we address a philosophical question raised by the model: if the 
predictions of our model are verified, should we then conclude that there is really a 
process of accumulation? Or is this merely a mathematical fiction? 

Physics has historically been built on analogies which have been proposed, explored, 
and pushed to their limits. For example. Maxwells equations were originally motivated 
by an analogy between electromagnetic fields and local displacements within an 
incompressible fluid medium due to stresses and strains acting upon that medium [8]. 
But as electromagnetic theory developed, the limitations of this analogy became 
increasingly apparentto the extent that it is scarcely mentioned in university courses 
on electromagnetism, and only a few vestiges may be seen in some of the terminology 
(such as stress tensor). Another important analogy (that has captured the popular 
imagination) is the idea that gravity bends space. This foundational idea motivated 
Einstein to look to differential geometry for mathematical formulations of the theory. 
The inability of general relativity to deal with quantum mechanics shows that the 
analogy can only go so far. More germane to the subject of this paper, the analogy 
between the statistical-mechanical partition function and the expression ([^ was a prime 
motivator for Euclidean quantum field theory [S]. 

The model in this paper is in some sense merely an extended analogy between 
quantum physics and wireless communications. It is our hope that further exploration 
of this analogy may lead to additional physical insights. Along these lines, we mention 
briefly some possibilities suggested by the model we have presented: 

• Feynman path integrals are well-known for yielding important physical results, 
while at the same time lacking mathematical rigor. Our model suggests making 
a correspondence between paths and possible states of the universe. This line 
of attack could lead to a less problematic mathematical characterization of path 
integrals. 

• Although formula ([^ is based on an action, so far we have said nothing about the 
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action, nor the fields that determine its value. Our analogy with signal processing 
suggests that there may be a relationship between the various types of quantum 
fields and signal modulations. 

• The model is designed to give an account of observed probability distributions for 
quantum events. However, so far we have not really defined event. Certainly this 
has something to do with the configuration of the fields involved: and perhaps this 
also may be understandable in terms of a signal-based representation of the fields. 

Although we hope that our model will be a source of insight, we recognize that even in 
the best case, the analogy that we have suggested will have proscribed limits and should 
not be thought of as representing “ultimate reality. Nonetheless, if this conceptual model 
proves to be accurate, it has profound implications for how we may regard the world 
around us, and how we regard ourselves as free agents within it. 


6. Simulation Code 


The following Octave/Matlab code was used for the simulation in Figures]^ and 


% Parameters 
clear all; 
nsim = 10000000; 
nconfig=l; 

Theta_fac = 10; 

Theta_fac0=Theta_fac; 
Ncfg = 31; 
n_acc_mean = 625; 
max_jump = 2; 
acc_mean0 = n_acc_mean; 


°/o # simulations 

70 # configs simulated 
% Theta increment 

"/o Orig. theta 

71 Number of internal configs 
°/o M interval 

7o For Markov — max jump 
7o Orig. M 


p = l/n_acc_mean; 

Theta = Theta_f ac*sqrt (n_acc_mean); 7o Rescaled threshold 


7o Arrays to store results 
Counts = zeros(Ncfg,1); 

Q = □; 


7o Create measurable configurations 
Psi = cos((0:1:Ncfg-1)/(Ncfg)*2*pi)'; 
Prob = abs(Psi)."2; 

Prob = Prob / sum(Prob); 


7o Computations 
for jj = l:l:nconfig 
for ii = l:l:nsim 


7o Loop over configurations 
7o Perform simulations 
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A = 0; 

this_cfg = randi([0,Ncfg-l]); % Choose current w 

°/o Accumulate: 

while abs(A) < Theta X Until threshold is attained 

this_cfg = mod(this_cfg + sign(randn)*randi([1,max_jump]),Ncfg); 
Ptmp = Psi(this_cfg+1);% Amplitude 
while randO > p 

Rtmp = randn()+li*randn(); 

A = A + Rtmp*Ptmp; 

if abs(A)>Theta "Z, If pass threshold, the break and record w 
break 

end 

end 

end 

Counts(this_cfg+l) = Counts(this_cfg+l)+l; % record w 

end 

Q = [Q Counts/sum(Counts)] % Summary results for this config 
Theta_fac = Theta_fac + Theta_facO; "Zlncrement theta 
Theta = Theta_fac*sqrt(n_acc_mean); 
end 

Prob = abs(Psi)."2; 

Prob = Prob / sum (Prob); "/oNormalized, sorted probabilities (for theory) 
plot(Prob,Q - Prob*ones(1,nconfig),'*O; 
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